Journal  of  Power  Sources  267  (2014)  566-575 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Journal  of  Power  Sources 

journal  homepage:  www.elsevier.com/locate/jpowsour 


Control  design  for  robust  tracking  and  smooth  transition  in  power 
systems  with  battery/supercapacitor  hybrid  energy  storage  devices^ 

Hoeguk  Jung,  Haifeng  Wang,  Tingshu  Hu 

Department  of  Electrical  and  Computer  Engineering,  University  of  Massachusetts,  Lowell,  MA  01854,  United  States 


CrossMark 


HIGHLIGHTS 


•  Control  design  problems  addressed  for  battery/supercapacitor  hybrid  energy  systems. 

•  Active  current  control  implemented  via  two  bi-directional  buck-boost  converters. 

•  Controller  designed  via  efficient  numerical  optimization. 

•  Robust  tracking  achieved  for  regulating  battery  current  and  dc-bus  voltage. 

•  Smooth  transition  achieved  during  load  switch. 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  7  February  2014 
Received  in  revised  form 
27  March  2014 
Accepted  7  May  2014 
Available  online  2  June  2014 


Keywords: 

Battery 

Supercapacitor 
Hybrid  energy  storage 
dc— dc  Converter 
Robust  tracking 
Smooth  transition 


This  paper  considers  some  control  design  problems  in  a  power  system  driven  by  battery/supercapacitor 
hybrid  energy  storage  devices.  The  currents  in  the  battery  and  the  supercapacitor  are  actively  controlled 
by  two  bidirectional  buck-boost  converters.  Two  control  objectives  are  addressed  in  this  paper:  one  is  to 
achieve  robust  tracking  of  two  reference  variables,  the  battery  current  and  the  load  voltage,  the  other  is 
to  achieve  smooth  transition  of  these  variables  during  load  switch.  Based  on  the  state-space  averaged 
model  we  newly  developed,  the  control  design  problems  are  converted  into  numerically  efficient  opti¬ 
mization  problems  with  linear  matrix  inequality  (LMI)  constraints.  An  experimental  system  is  con¬ 
structed  to  validate  the  control  design  methods. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Batteries  and  supercapacitors  have  been  combined  into  hybrid 
energy  storage  devices  which  have  both  the  high  energy  density  of 
the  batteries  and  the  high  power  density  of  the  supercapacitors. 
The  main  function  of  supercapacitors  in  a  hybrid  energy  storage 
system  is  to  provide  high  currents  during  hard  transience  such  as 
motor  start,  which  is  essential  to  protecting  the  batteries  from  fatal 
damages  caused  by  over-discharge.  Supercapacitors  are  also  used  to 
absorb  excess  current  when  wind/solar  power  is  abundant,  and  to 
store  energy  from  regenerative  braking.  The  battery/supercapacitor 
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hybrid  energy  storage  systems  have  been  widely  used  in  electric, 
hybrid  and  plug-in  hybrid  electric  vehicles,  e.g.,  see  Refs.  [1—4], 
They  have  also  found  applications  in  wind  systems  [5—7], 
communication  systems  [8],  photovoltaic  systems  [9],  and  micro¬ 
grids  [10].  In  some  hybrid  electric  vehicle  applications,  the  hybrid 
energy  storage  system  also  include  fuel  cells  as  one  power  source 
[11—13], 

There  are  various  configurations  to  combine  batteries  and 
supercapacitors.  The  advantage  and  disadvantage  of  each  configu¬ 
ration  are  discussed  in  Refs.  [3,14—16],  In  Ref.  [17],  the  performance 
of  different  configurations  were  compared  via  simulation  on  a 
stand  along  power  system.  The  simplest  configuration  is  to  connect 
the  battery  and  the  supercapacitor  in  parallel.  This  will  reduce  the 
current  stress  of  the  battery  but  the  two  power  sources  need  the 
same  voltage  and  the  power  flow  cannot  be  controlled.  To  achieve 
active  control  of  currents,  especially  the  current  from  the  battery, 
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dc— dc  converters  are  needed  to  connect  the  power  sources  and  the 
load  or  dc  bus.  In  some  configurations  [18—21],  one  dc— dc  con¬ 
verter  is  used  to  connect  the  battery  and  the  supercapacitor. 

The  most  commonly  used  configuration  contains  two  bi¬ 
directional  buck— boost  converter  connected  in  parallel  at  the 
load  side,  and  fed  by  a  battery  and  a  supercapacitor  respectively,  see 
Fig.  1.  The  advantage  of  this  configuration  is  that  both  the  current 
from  the  battery  and  the  current  from  the  supercapacitor  can  be 
actively  controlled.  Such  a  configuration  has  been  considered,  for 
example,  in  Refs.  [8,14,16,22].  Similar  configurations  have  been 
used  to  combine  fuel  cells  and  supercapacitors,  e.g.,  in  Ref.  [11], 
where  the  fuel  cell  is  connected  to  the  dc  bus  via  a  one-directional 
boost  converter.  In  Ref.  [23],  three  buck— boost  converters  were 
connected  in  parallel  to  combine  three  power  sources.  For  the 
configuration  in  Fig.  1,  different  strategies  have  been  proposed  to 
actively  control  the  current  flow  from  the  battery  and  the 
supercapacitor. 

A  common  strategy  is  to  use  a  certain  power  management  al¬ 
gorithm  to  determine  a  reference  current  that  is  needed  from  the 
battery  or  the  supercapacitor,  then  use  decentralized  reference 
tracking  control,  typically  PI  control,  on  each  dc— dc  converter  to 
track  the  respective  reference  current.  In  Ref.  [24  ,  one  algorithm 
was  proposed  for  maximum  efficiency  and  the  other  algorithm  for 
minimal  battery  current  to  prolong  the  life  span  of  batteries.  In 
Ref.  [16],  a  power  management  algorithm  was  proposed  for  mini¬ 
mization  of  the  magnitude/fluctuation  of  the  battery  current  and 
energy  loss.  In  Ref.  [25],  three  piecewise  linear  functions  were 
proposed  for  allocating  the  battery  current  and  supercapacitor 
current.  Simulation  results  were  generated  to  demonstrate  the 
effectiveness  of  the  allocation  function  under  different  operating 
conditions.  In  Ref.  [10],  several  current  allocation  strategies  were 
proposed  by  considering  the  power  demand,  the  state  of  charge  of 
the  battery  and  the  supercapacitor.  In  Ref.  [22],  three  classical 
control  loops  were  constructed  to  regulate  the  battery  current,  the 
output  voltage  and  the  supercapacitor  voltage. 

In  most  of  the  existing  literature  which  considered  similar 
configuration  as  in  Fig.  1,  the  control  design  algorithms  were 
focused  on  generating  reference  currents  for  the  battery  and  the 
supercapacitor.  It  was  generally  assumed  that  a  classical  reference 
tracking  loop  would  yield  a  required  current  from  the  battery  or  the 
supercapacitor.  These  control  strategies  do  not  consider  the  inter¬ 
action  among  the  control  loops  and  the  power  loss  in  the  circuit 
elements.  In  fact,  the  parallel  dc— dc  converters  are  coupled  via  the 


same  dc-bus  or  load  and  have  complex  interactions.  If  the  in¬ 
teractions  among  the  different  control  loops  are  ignored,  the 
transient  behavior  may  not  be  desirable,  such  as  with  large  over¬ 
shoot  or  undershoot,  which  may  damage  the  circuit  parts.  If  the 
power  loss  is  not  carefully  considered  in  the  current  allocation  al¬ 
gorithm,  the  output  power  may  be  different  than  the  desired 
power. 

For  better  analysis  and  control  of  the  coupled  dc— dc  converters, 
we  need  effective  and  faithful  models  for  the  whole  power  elec¬ 
tronic  system.  Recently,  we  derived  two  state-space  averaged 
models  for  the  system  in  Ref.  [26],  One  for  simulation  and  analysis 
of  the  open-loop  system  and  the  closed-loop  system,  and  the  other 
for  control  design.  The  advantage  of  the  averaged  models  is  the 
simplicity  and  that  they  have  clear  structure  in  terms  of  the  two 
duty  cycles.  Furthermore,  they  are  much  faster  for  simulation  than 
the  SimPower  models. 

In  this  paper,  we  would  like  to  address  some  control  design 
problems  for  the  typical  configuration  with  two  bidirectional 
dc— dc  converters  as  in  Fig.  1.  A  standard  control  problem  in  such 
power  systems  is  to  achieve  reference  tracking  for  some  variables 
such  as  battery  current,  supercapacitor  current,  load  voltage  and 
load  current.  Since  the  hybrid  energy  storage  system  has  two 
control  inputs,  the  duty  cycles  of  the  two  dc— dc  converters,  it 
should  be  able  to  track  references  for  two  circuit  variables.  One  may 
still  choose  the  battery  current  and  the  supercapacitor  current  as  in 
some  earlier  works.  In  this  paper,  we  will  choose  the  battery  cur¬ 
rent  and  the  load  voltage  since  the  supercapacitor  is  playing  a 
supporting  role  and  it  can  supply  or  absorb  almost  any  current  as 
needed.  After  all,  the  state-of-health  of  the  battery  (which  depends 
on  its  charging/discharging  profiles)  and  the  performance  of  the 
load  are  of  high  priority. 

We  will  also  address  the  problem  of  control  design  for  smooth 
transition  of  some  key  variables  during  load  switch.  This  will  also 
help  to  protect  the  battery  and  to  enhance  the  performance  of  the 
whole  system. 

It  should  be  noted  that  this  paper's  control  design  methods  can 
be  combined  with  the  existing  current  allocation  strategies.  With 
the  reference  currents  generated  by  any  of  the  existing  method,  the 
tracking  control  can  be  implemented  by  the  controllers  designed  by 
our  optimization  algorithm. 

This  paper  is  organized  as  follows.  In  Section  2,  we  describe  the 
circuit  system  and  briefly  summarize  the  two  state-space  averaged 
models  developed  in  Ref.  [26].  In  Section  3,  two  design  objectives 


Fig.  1.  Parallel  topology  of  the  buck-boost  converters. 
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are  addressed:  One  is  to  achieve  robust  tracking  and  the  other  is  to 
achieve  robust  tracking  with  smooth  transition  during  load  switch. 
Both  design  problems  are  formulated  as  optimization  problems 
with  linear  matrix  inequality  constraints.  In  Section  4,  the  design 
methods  are  applied  to  an  experimental  system  and  yield  two  state 
feedback  laws.  The  performances  of  the  system  under  these  feed¬ 
back  laws  are  validated  via  the  averaged  model  and  the  SimPower 
Model.  Section  5  validates  the  effectiveness  of  the  control  design 
methods  on  an  experimental  system.  Section  6  concludes  this 
paper. 


2.  Description  of  the  power  system  and  two  state-space 
averaged  models 

The  circuit  diagram  of  the  hybrid  energy  storage  system  is 
illustrated  in  Fig.  1.  It  consists  of  two  standard  bi-directional 
buck— boost  converter  connected  in  parallel  at  the  load  side,  and 
fed  by  a  battery  and  a  supercapacitor  respectively.  The  load  is 
represented  by  a  resistor  for  simplicity,  as  in  some  literature,  e.g., 
[27],  In  Ref.  [22],  the  load  is  represented  by  a  current  source.  In 
electrical  vehicle  applications,  the  load  is  typically  an  inverter  fol¬ 
lowed  by  a  motor  [18],  The  inverter  and  the  motor  can  be  made 
equivalent  to  a  resistor  whose  resistance  varies  within  a  certain 
range.  During  motor  start,  the  equivalent  load  resistance  is  much 
smaller  than  in  other  operating  conditions. 

For  simplicity  of  presentation,  the  battery  and  the  super¬ 
capacitor  are  described  with  models  which  only  include  parasitic 
series  resistances  Re  and  Ru,  respectively.  More  comprehensive 
models  can  be  considered,  such  as  adding  one  or  two  pairs  of 
parallel  resistor  and  capacitor  for  the  battery,  as  in  Refs.  [28,29],  and 
using  a  second-order  RC  circuit  to  describe  a  supercapacitor,  as  in 
Refs.  [1,27],  In  Ref.  [28],  an  algebraic  method  was  developed  for 
identifying  the  circuit  parameters  of  a  battery. 

The  two  MOSFETs  Si  and  S 2  operate  synchronously,  when  one  is 
off,  the  other  is  on.  Similarly  with  S3  and  S4.  Let  the  duty  cycle  for 
the  on  state  of  Si  be  D 1  and  the  duty  cycle  for  the  on-state  of  S3  be 
D2.  The  on-resistances  for  Si,  S2,  S3,  S4  are  denoted  as  R0n,i,  R0n,2. 
Ron, 3.  Ron, 4,  respectively.  Based  on  the  state-space  averaging 
method,  we  derived  two  state-space  averaged  models  in  Ref.  [26] 
for  the  hybrid  energy  storage  system:  a  6th  order  model  for 
simulation  and  a  5th  order  model  for  control  design.  If  more 
complex  models  for  the  battery  and  the  supercapacitor  are 
considered,  the  method  in  Ref.  [26]  will  yield  higher  order  state- 
space  models,  with  matrices  of  the  same  structure  but  larger 
dimensions. 

For  the  6th  order  model,  the  voltage  vu  of  the  supercapacitor  is 
considered  as  a  state  variable  in  addition  to  Vi,  V2,  h,  i2  and  v0,  as 
marked  in  Fig.  1.  The  6  state  variables  are  arranged  in  a  state  vector 
£  =  [tq  v2  vu  ii  i2  v0]T.  The  averages  of  the  state  variables 
over  one  switching  period  are  denoted  as  v\ .  V2,  vu,  h-  h,  F0, 
respectively.  The  average  of  the  state  vector  £  is  denoted  as  C. 

The  6th  order  averaged  model  is 

?  =  (Ao  +  A^t  +  A2D2^£  +  BeE  (1) 

where  E  is  the  open  circuit  voltage  of  the  battery.  The  matrices 
A0,  Ai  ,  A2  and  Be  depend  on  circuit  parameters.  Detailed  structure 
and  expression  for  these  matrices  can  be  found  in  Ref.  [26], 

It  is  commented  in  Ref.  [26]  that  the  6th  order  model  (1)  is 
suitable  for  simulation  of  the  open-loop  system  under  constant 
duty  cycles  Di  and  D2,  as  well  as  for  the  closed-loop  system  under 
particular  control  strategies.  However,  it  is  not  suitable  for  steady- 
state  analysis  or  control  design.  The  main  reason  is  that,  its 
equilibrium  point,  where  £  =  0,  is  not  a  useful  operating  condition 


since  it  implies  Avu/A t  —  0.  Under  this  condition,  the  super¬ 
capacitor  is  not  supplying/absorbing  current  and  thus  not  assisting 
the  power  system.  In  general  application,  the  supercapacitor  is 
used  for  discharging  high  current  for  only  a  brief  moment  and  will 
be  recharged  soon.  The  high  current  discharging  mode  is  very 
important  but  will  not  last  long  enough  to  reach  a  steady  state. 

For  control  design  (stabilization  or  tracking),  we  need  to 
consider  a  nominal  working  condition  and  to  derive  a  perturba¬ 
tion  model  around  this  nominal  condition,  which  has  to  be  a 
steady  state  (or  an  equilibrium  point).  To  handle  this  situation,  we 
may  replace  the  supercapacitor  Cu  with  an  ideal  voltage  source  Eu 
for  the  purpose  of  control  design  and  for  determining  the  input/ 
output  current  or  voltage  after  the  fast  transience  of  the  dc— dc 
converters.  We  can  ignore  the  dynamics  of  the  supercapacitor 
even  when  it  is  discharging  at  a  high  current,  since  the  decrease  of 
its  voltage  is  still  slow  as  compared  to  the  fast  dynamics  of  other 
elements. 

In  this  respect,  a  5th  order  averaged  model  were  derived  in 
Ref.  [26].  Denote  the  state  variable  as  %=  [v\  v2  q  i2  v0]T,  and 
the  averaged  state  as  f.  Denote 
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Denote  Ve  = 


,  then  the  5th  order  averaged  model  is 


?  =  (A)  +A}Di  +  A2D2^ 


+  BeVe 


(5) 


3.  Control  design  for  robust  reference  tracking  and  smooth 
transition 


In  this  paper  we  focus  on  the  control  of  the  fast  dynamics  charac¬ 
terizing  the  power  converters,  neglecting  the  above-mentioned 
supercapacitor  sizing  and  voltage  regulation.  Slow  variation  of 
the  supercapacitor  voltage  will  be  regarded  as  a  model  uncertainty 
to  be  handled  by  the  robust  feedback  controller. 

Combining  the  5th  order  model  (5)  with  an  output  of  interest, 
we  obtain 


3.1.  Converting  the  tracking  problem  to  a  stabilization  problem 


Strictly  speaking,  stabilization  for  the  original  hybrid  energy 
storage  system  as  described  Fig.  1  and  the  6th  order  averaged  model 
( 1 )  is  not  a  meaningful  problem.  This  is  because  stability  or  tracking 
is  defined  in  terms  of  a  useful  equilibrium  point,  where  you  want  the 
state  of  the  system  to  stay  there  for  a  desirable  time  period.  The 
equilibrium  point  of  the  6th  order  model  is  useless  since  the  de¬ 
rivative  of  the  supercapacitor  voltage  must  be  0,  which  means  that 
its  current  is  0.  However,  we  need  a  perturbation  model  around  the 
equilibrium  point  for  control  design.  So  we  have  to  use  the  5th  order 
model,  which  is  obtained  by  excluding  the  voltage  of  the  super¬ 
capacitor  from  the  state  of  the  6th  order  model,  for  control  design. 
Instead,  the  supercapacitor  voltage  is  considered  as  a  voltage  source 
Eu,  which  is  allowed  to  be  varying  "slowly"  within  a  certain  range. 

The  controller  designed  for  the  5th  order  model  should  work  for 
the  6th  order  model  for  reasonably  long  time  period.  The  reasons 
include:  1)  the  controller  designed  by  our  LMl-based  method  is 
robust,  which  means  that  it  can  tolerate  some  uncertainty  or 
change  of  parameters;  2)  The  change  of  supercapacitor  voltage  is 
very  slow  as  compared  to  the  transient  response  of  other  state 
variables.  The  sizing  of  the  supercapacitor  should  take  into  account 
this  requirement. 

Since  there  are  two  control  inputs  D 1  and  D 2,  a  practical  control 
objective  is  to  track  the  reference  for  a  two-dimensional  output, 


e-g-,y  = 


V 

h 

.'sc. 

Vo  . 

.  In  general,  consider  y  =  Cf  where  Cis  a  matrix 


of  two  rows. 

For  safety  and  efficiency,  the  duty  cycles  Di,  D 2  are  restricted  in  a 
subset  of  (0,  1).  For  this  reason  and  other  circuit  limitations,  we 
make  the  following  observations: 


1.  For  the  5th  order  model  where  the  supercapacitor  is  replaced 
with  an  ideal  voltage  source  Eu,  for  each  output  reference  yref, 
there  is  a  certain  range  for  the  source  voltages  (£,  Eu)  where 
tracking  is  possible. 

2.  According  to  charge  and  energy  conservation  principles,  for  the 
original  6th  order  system  with  supercapacitor,  generic  tracking 
can  only  last  for  a  finite  time  period,  beyond  which  the  super¬ 
capacitor  voltage  will  drop  (or  rise)  out  of  the  range  where  yref 
can  be  tracked. 


To  see  item  1,  we  note  that  v0  has  to  be  greater  than  both  E  and  Eu 
and  the  gain  from  E  and  Eu  to  v0  is  limited  by  the  load  resistance  and 
the  range  of  duty  cycles.  For  item  2,  we  note  that  in  a  real  hybrid 
energy  storage  system,  the  reference  yref  is  recalculated  after  a  while 
according  the  load  requirement,  the  state  of  charge  of  the  battery  and 
the  state  of  charge  of  the  supercapacitor.  During  heavy  load  period, 
the  supercapacitor  supplies  most  of  the  current  and  its  voltage  de¬ 
creases  quickly.  This  period  does  not  (and  cannot)  last  long.  During 
light  load  period,  the  supercapacitor  is  recharged.  The  charging  rate 
depends  on  the  battery  current.  This  period  can  last  much  longer  than 
the  heavy  load  period.  But  if  the  same  operation  mode  is  allowed  to 
last  forever,  the  supercapacitor  voltage  will  rise  out  of  the  limit. 

Item  2  in  the  above  statements  further  justifies  the  use  of  the 
5th  order  model  for  control  design  instead  of  the  6th  order  model. 


f  —  (A)  +A1D1  +  A2D2)  ?  +  BeVe 
y  =  c? 


where  Ve  = 


and  C  is  a  2  by  5  matrix  depending  on  the  circuit 


variables  to  be  tracked.  For  example, 
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For  a  particular  pair  (£,  £u),  a  reference  yref  can  be  tracked  if 
there  exist  Dj,  D2e( 0.  1)  such  that  Cfss  =  yref  where  Jss  is  the 
steady-state  satisfying  |ss  =  0,  i.e., 
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Before  designing  a  control  law,  we  need  to  choose  a  suitable 
nominal  operating  condition  which  involves  the  following  vari- 
'£0 


ables:  Ve0  = 


cu  0 


,  D10,  D20,  and  ?ss,0  satisfying 


?ss,0  =  —  (A)  +  AD10  +  A2D2(j) 


BeVeo 


(8) 


Next,  we  need  a  perturbation  model  around  this  nominal  con¬ 
dition.  Define  x  =  f  -  fss,0,  y  =  y  -  C?ss0,  10  =  Dx  -  Dw, 


u2  —  D2  —  D2 0  and  u  = 


.  Also  denote 


A)  =A)  +a1d10  +a2d20,  b=  [a?SSi0  A2?SSio] 


then  by  plugging  ?  =  x  +  fSs,o.  Di  =  Ui +  D10  and  D2  =  u2  +  D20  into 
(6)  and  applying  (8),  we  obtain  the  following  perturbation  model: 

x  =  A)X  +  A]XUi  +A2xu2  +Bu  +  Be(Ve  -  Ve0)  rq-, 

y  =  Cx  K  1 

We  may  design  a  feedback  law  to  stabilize  the  system  (9)  under 
the  nominal  condition  where  Ve  =  Veo.  However,  it  is  almost  always 
true  that  Ve  =£  Veo  since  the  battery  voltage  or  the  supercapacitor 
voltage  are  always  changing.  Furthermore,  the  desired  value  for  the 
output  is  also  changed  frequently  and  the  load  resistance  is  not  a 
constant. 

A  common  strategy  to  achieve  robust  reference  tracking  in  the 
presence  of  uncertainties,  is  to  integrate  the  tracking  error 
e=y  -  yref  and  include  the  integration  as  a  new  state  variable.  Here 
yref  is  a  two  dimensional  vector.  It  is  generally  different  than  C?ss  0, 
the  output  at  the  nominal  condition.  To  be  specific,  define 


(10) 


and  define  the  augmented  state  as  xw  := 
following  matrices: 


x 

xa 


Also  define  the 
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A  = 
B  = 


A0  0 
C  0 
B 


>  Aftl  = 


0 

0  0 


i  A/,2 


0 


,  s 


Be(ve  -  Ve0) 
-yref 


A2  0 

o  o 


where  the  zero  blocks  have  compatible  dimensions.  Then  xw 
satisfies 


Xw  —  AXw  +  AfoiXwU]  H-  Afo2XwU2  +  Bu  -j-  g  (1 1 ) 

Here  xweR7  includes  the  perturbations  of  the  voltages  of  three 
capacitors  Cj,  C2,  C0,  the  two  inductor  currents,  and  two  integrator 
outputs.  This  turns  the  reference  tracking  problem  into  a  standard 
stabilization  problem:  design  a  feedback  control  law  u  =  f[xw) 
which  stabilizes  (11)  with  a  large  stability  region,  under  the  nom¬ 
inal  condition  g=  0,  and  under  the  control  constraint  D\,  D2e(0, 1). 
If  the  system  drifts  away  from  the  nominal  working  condition  and 
goes  to  another  equilibrium  point,  each  state  variable,  in  particular, 
the  integral  xa,  will  still  reach  a  steady  state.  This  means  thaty  -  yref 
must  go  to  zero  and  the  output  y  is  regulated  to  the  desired  value 
yref- 

In  what  follows,  we  design  a  stabilizing  feedback  for  (11)  under 
input  constraint  by  extending  the  method  in  Refs.  [29,30]  for  sys¬ 
tems  with  one  input  to  systems  with  two  inputs. 


3.2.  State  feedback  design  for  robust  tracking  via  LMI  optimization 


Let  us  first  deal  with  the  input  constraints.  In  general,  assume 
Dje[Djmin.  Djmax]  c[0,  1],  j  =  1,  2.  This  corresponds  to 
H/mjn  —  Dp  <  Uj  <  Djmax  —  Djq.  Denoting  Umj  —  D]0  —  Djm  jn, 
uPj  =  Djmax  -  Djo,  the  input  constraints  can  be  expressed  as 

—Umj  —  Uj  —  t/pj,  ./  -  1  -  2  (12) 

The  input  constraint  will  be  enforced  via  some  saturation 
functions  defined  as  follows.  For  j  =  1,  2 

f  UPJ  uJ>uPj 

Sat(lfj)  =  Uj  Uj  e  [  -  umj,  upj] 

[  -Umj  Uj  <  -  Umj 


with  some  abuse  of  notation,  for  ueR2,  denote  sat(u)  = 


sat(Uj) 
sat(u2) 

We  will  consider  a  simple  saturated  state  feedback  u  =  sat(Xx). 
Under  this  state  feedback  and  the  nominal  working  condition,  we 
obtain  the  following  closed-loop  system: 


feedback  gain  I<  to  achieve  the  maximal  stability  region  under  the 
input  constraint  (12).  Note  that  the  size  of  the  stability  region  is 
closely  related  to  the  robustness  against  parameter  change  and  the 
range  of  references  that  can  be  tracked. 

The  design  methods  are  developed  based  on  Lyapunov  function 
and  some  theory  on  set  invariance,  which  are  powerful  tools  for 
dealing  with  input  and  state  constraints.  The  stabilization  problem 
can  be  converted  into  a  numerically  efficient  optimization  problem 
with  linear  matrix  inequality  (LMI)  constraints.  Detailed  ideas  and 
explanations  can  be  found  in  Refs.  [30—32]. 

For  the  particular  system  (14)  under  input  constraint  (12),  the 
design  objective  for  achieving  the  maximal  stability  region  can  be 
cast  into  the  following  optimization  problem: 


inf  y 

q>o,y,t 

s.t.  a)AjQ_  +  qAJ  +  BY+  YTBT  < 


2vQ.,  i  =  1,  ...,4 


b) 


c) 


mm(u2mj.u2pj 


YJ 


Q  I 

I  y  I 


>  0 


Yi 

Q 


>0,  j  =  l,2 


(15) 


where  Yf,  Y2  are  the  two  rows  of  Y  and  77  >  0  is  an  optimization 
parameter  that  can  be  chosen  by  experience.  In  general,  larger  77 
would  result  in  larger  elements  in  the  feedback  gain  K  and  faster 
response.  However,  if  the  elements  in  K  are  too  large,  the  control 
input  will  be  easily  saturated,  which  may  yields  small  stability  re¬ 
gion.  Thus  we  usually  have  to  adjust  77  after  some  simulation  or 
experiment  on  the  closed-loop  system. 

The  matrix  variables  to  be  optimized  are  Q  eR7x7  and  Y eR2x7, 
while  the  other  matrices  are  constants  obtained  from  the  averaged 
model  following  some  steps.  The  optimization  problem  can  be 
easily  solved  with  the  “mincx”  solver  (minimizing  linear  objective 
under  LMI  constraints)  in  the  LMI  toolbox  of  MATLAB.  After  Qand  Y 
are  solved,  the  feedback  matrix  is  obtained  as  I<  =  YQ 

The  feedback  control  designed  via  the  LMI  optimization  prob¬ 
lem  will  always  allow  certain  ranges  of  parameters  where  robust 
tracking  can  be  achieved.  These  ranges  can  also  be  estimated  via 
some  optimization  problem  but  the  results  may  be  conservative. 
The  averaged  model  will  provide  an  easy  and  fast  approach  for 
verifying  robustness  within  a  certain  parameter  range. 


3.3.  State  feedback  design  for  smooth  transition 


xw=  (a  +  Ablsat(/C]Xw)  +  Ab2sat(k2xw)^xw  +  Bsat(/<xw)  (13) 

The  nonlinear  matrices  A^satjfcjXw)  and  Ab2sat(k2Xw)  will  be 
handled  with  inclusion  and  the  term  Bsat(/<xw)  will  be  handled  by 
the  method  in  Ref.  [31],  Define 

Al  =  A  -  umlAbl  -  um2Ab2.  A2  —  A  —  umiAbl  +  up2Ab2 
A3  =  A  +  up]Abl  -  um2Ab 2,  A4  =  A  +  up2Abl  +  up2Ab2 

Then  xw  satisfies 

xw e co|a,xw  +  Bsat(/<xw) ,  i  =  l,  ...,  4j  (14) 

where  co{X}  denote  the  convex  hull  of  the  set  X.  In  particular, 
co{a,  b.c,  d}  =  {y1a  +  y2b  +  y3C  +  y4d  :  =  1,  yt  >  0} 

By  applying  the  methods  reported  in  Refs.  [30— 32  ,  to  the  dif¬ 
ferential  inclusion  (14),  with  slight  generalization,  we  can  design  a 


Simulation  and  initial  experimental  results  showed  that  the 
controller  designed  for  robust  reference  tracking  may  yield  unsat¬ 
isfactory  transient  performances.  For  example,  when  the  load  is 
switched  from  200  Q  to  4  fi,  there  is  significant  undershoot  for  the 
load  voltage  and  large  overshoot  for  the  battery  current  (see  Fig.  2). 
When  the  load  is  switched  from  heavy  to  light,  there  are  significant 
overshoot  for  load  voltage  and  undershoot  for  battery  current. 
These  overshoots  and  undershoots  will  have  adverse  effects  to  the 
battery,  the  load  performance  and  the  reliability  of  the  whole 
system. 

To  reduce  the  overshoot/undershoot  of  the  outputs,  we  first 
need  some  tools  to  evaluate  such  transient  performances.  Such 
tools  were  developed  in  Ref.  [33  ,  where  the  overshoot/undershoot 
is  evaluated  via  invariant  ellipsoids.  Invariant  ellipsoid  is  a  powerful 
tool  based  on  Lyapunov  theory.  If  there  is  an  invariant  set  S  that 
contains  all  the  possible  initial  conditions,  then  the  state  trajec¬ 
tories  will  stay  within  the  invariant  set  and  we  can  estimate  the 
maximal/minimal  output  via  the  invariant  set.  Suppose  that  the  set 
of  initial  conditions  is  bounded  by  a  ball  Xq  =  {xeRn  :  xTx  <  r^}, 
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Response  under  K  ,  E=12V,  v  (0)=14V,  v  =20V,  i  =1.5A 

r  1  u'  '  o.ref  b.ref 


time  (second) 


Fig.  2.  Response  to  load  switch:  £  =  12  V,  v„(0)  =  14  V,  feedback  gain 

and  the  invariant  ellipsoid  is  e  =  {x£R“  :  xT()_1x  <  1}.  The  rela¬ 
tionship  X0  ce  can  be  described  by  the  matrix  inequality  Q  >  rpn. 

Suppose  that  we  need  to  evaluate  the  maximal  y \  =  C\x  under 
the  condition  x0eX0.  Assume  X0ce  and  e  is  an  invariant  set.  A 
bound  on  the  maximal  yi  can  be  computed  as 


inf 7^,  s.t. 


(16) 


To  design  a  state  feedback  law  with  minimal  overshoot  for  a 
certain  output,  we  just  need  to  impose  the  LMI  constraint  in  (16) 
together  with  the  constraints  in  (15),  with  the  objective  of  mini¬ 
mizing  yi  under  a  guaranteed  y. 

For  the  particular  system  in  this  paper,  we  need  to  reduce  the 
overshoot/undershoot  of  two  outputs:  yi  =  Cjx  as  the  battery 
current  andy2  =  C2X  as  the  load  voltage.  We  use  yi  and  72  to  denote 
the  magnitudes  of  yi  and  y2  respectively. 

Since  these  two  outputs  may  have  conflict,  we  can  use  a 
weighting  coefficient  w  and  consider  a  combined  measure 
y\  +  wy\.  The  optimization  problem  can  be  described  as 


inf  y\  +  wy  2 


-T  - 


Q>  0X71,72 

s.t.  a)AiQ  +  QA‘i  +  BY  +Y‘B‘  < -2pQ_.  i  =  1 , . . . ,  4 


■=r 


b) 

c) 

d ) 


«■  -upj) 

Yj 

Y? 

Q 

>0,J  =  1,2 


(17) 


Q  I 

I  yi 


QC, 


>  0 


CiQ 

>  0 

72 

C2Q" 

Q 

qc2t 

Q 

>  0. 


As  compared  with  the  optimization  problem  (15),  the  new 
constraint  d)  is  added  with  an  objective  of  minimizing  y |  +  wy\.  To 
ensure  a  certain  robustness  on  the  tracking  performance,  we 
choose  an  appropriate  value  for  y. 

Note  that  there  are  three  parameters,  r],  y  and  w,  that  need  to  be 
chosen  before  optimization.  Since  the  relationship  between  these 
parameters  and  the  output  performances  are  complex,  we  may  try 
different  combinations  and  then  use  simulation  to  decide  which 
resulting  controller  is  the  best. 


4.  Control  design  for  an  experimental  system  and  simulation 
results 


4.3.  Open-loop  system  description 


An  experimental  system  was  constructed  for  the  hybrid  energy 
storage  configuration  as  described  in  Fig.  1.  The  fixed  parameters 
are  provided  in  Table  1. 

The  load  resistance  R  is  variable.  A  time  varying  profile  of  the 
resistance  R  can  be  constructed  to  mimic  the  different  operation 
modes  of  an  electric  machine.  In  our  simulation  and  experiment, 
we  used  a  simple  switch  to  change  the  resistance  between  a  few 
ohms  (heavy  load)  and  several  hundred  ohms  (light  load).  The 
battery  is  a  lead-acid  one  rated  12  V,  13  Ah.  Its  open  circuit  voltage 
varies  between  11.6  V  and  12.6  V,  corresponding  to  10%  and  100% 
state  of  charge,  respectively.  The  supercapacitor  consists  of  two 
parallel  ones  rated  58  F,  16.2  V  with  serial  resistance  0.022  0.  For 
control  design,  the  supercapacitor  was  considered  as  an  ideal 
voltage  source  varying  between  12  V  and  16  V.  It  should  be  clarified 
that  in  both  simulation  and  experiment,  the  supercapacitor  was 
used.  The  simulation  was  carried  out  with  both  SimPower  and 
Simulink  with  the  6th  order  averaged  model  (1 )  instead  of  the  5th 
order  model. 

The  control  objective  is  to  keep  the  load  voltage  at  a  desired 
value,  e.g„  20  V,  for  both  heavy  load  and  light  load  conditions.  To 
ensure  safe  discharging  of  the  battery,  we  would  like  to  keep  the 
battery  current  at  a  certain  constant  below  3  A.  The  reference  value 
for  the  battery  current  can  be  changed  according  to  the  state  of 
charge.  Thus  we  need  to  design  a  controller  for  tracking  the  output 


y  = 


>1 

Vo 


The  nominal  working  condition  is  chosen  based  on  the  control 
objective  and  the  rating  of  the  battery  and  the  supercapacitor.  The 
nominal  battery  voltage  is  chosen  as  E 0  =  12  V.  For  control  design, 
the  nominal  voltage  for  the  ideal  voltage  source  in  place  of  the 
supercapacitor  is  chosen  as  Eu 0  =  14  V.  For  the  load  resistance,  we 
focus  on  the  heavy  load  condition  thus  we  choose  R0  =  4  O.  The 


nominal  condition  for  the  output  is  chosen  as  y0 


3 

20 


.  The 


unique  duty  cycles  (Dh  D2 )  that  produce  this  nominal  output  is  (Dio, 
D20)  =  (0.447,  0.337)  and  the  steady  state  for  the  5th  order  averaged 
model  is  |ss0  =  [11.88  13.94  5.03  20]r.  Using  these  condi¬ 
tions,  a  perturbation  model  can  be  computed. 

For  the  duty  cycle,  we  impose  the  restriction  Di,  D2  e  [0,  0.8],  The 
lower  bound  0  and  upper  bound  0.8  are  used  to  compute 


timi  =  0.447,  Up  1  =  0.353,  um 2  =  0.337,  up 2  =  0.463. 


4.2.  Controller  for  robust  tracking  and  simulation 


First,  we  found  the  feedback  gain  I<  to  achieve  robust  tracking. 
By  solving  (15)  with  different  r\,  we  can  obtain  a  family  of  feedback 
gains  K.  Here  is  a  feedback  gain  corresponding  to  rj  -  50: 


-0.0005 

-0.0000 

0.0000 

-0.0000 

-0.0082 

-0.0001 

0.0000 

-0.0002 

0.0007 

-0.0009 

-2.8254 

0.6908 

-1.9620 

-2.8779 

This  feedback  gain  was  applied  to  both  the  averaged  model  and 
the  SimPower  model  (with  switching  frequency  20  kHz).  The  nominal 


reference  yref 


3 

20 


was  successfully  tracked  in  both  models. 
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Table  1 

Circuit  parameters. 


Symbol 

Description 

Quantity 

Unit 

£1 

Battery  side  inductance 

680 

[pH] 

Re 

Battery  series  resistance 

0.04 

[Q] 

Rn 

Inductor  series  resistance 

0.25 

[Q] 

Cl 

Battery  side  filter  capacitance 

1 

[mF] 

£.2 

Supercapacitor  side  inductance 

39 

[pH] 

R u 

Supercapacitor  serial  resistance 

0.011 

[Q] 

Rl2 

Inductor  series  resistance 

0.114 

[Q] 

C2 

Supercapacitor  side  filter  capacitance 

0.22 

[mF] 

Csc 

Supercapacitor  capacitance 

116 

[F] 

Ron 

MOSFET  on-resistance 

0.021 

[Q] 

c 

Load  side  capacitance 

1 

[mF] 

It  should  be  noted  that  the  simulation  with  the  averaged  model 
is  much  faster  than  that  by  the  SimPower  model.  For  the  averaged 
model,  the  sampling  time  can  be  chosen  as  10~4  s.  For  the  Sim¬ 
Power  model,  the  sampling  time  has  to  be  much  smaller  than  the 
switching  period.  Generally,  there  should  be  about  25  or  more 
samples  within  a  switching  period.  For  a  20  kHz  switching  fre¬ 
quency,  the  sampling  period  for  the  SimPower  model  needs  to  be 
2  x  10~6  s.  For  example,  to  run  a  simulation  of  0.6  s  as  in  Figs.  2—5, 
the  averaged  model  requires  18  s,  while  the  SimPower  model  re¬ 
quires  960  s.  Thus  the  averaged  model  can  save  a  lot  of  time  for 
determining  the  range  of  parameters  for  which  tracking  is 
achievable.  If  tracking  is  successful  for  a  particular  set  of  parameters 
by  the  averaged  model,  we  then  use  SimPower  model  to  validate  it. 

To  validate  the  robustness  of  the  tracking  performance,  we 
changed  the  load  resistance  R,  the  reference  values  for  the  output 
yref.  the  battery  voltage  E  and  the  initial  voltage  of  the  super¬ 
capacitor  vu(0)  and  run  the  simulation. 

Simulation  results  showed  that  robust  tracking  can  be  achieved 
within  a  large  range  of  parameters.  For  example,  with  R  =  4£l  and 


Tret  “ 


1.5 

20 


the  tracking  can  be  achieved 


for 


£e[ll,  19]  V,  vu(0)e[9,  19]  V.  For  larger  resistance  R,  larger 
ranges  for  £  and  vu(0)  are  allowed.  If  R  is  reduced  to  2  Q,  then  the 
source  voltages  cannot  be  too  low:  £  >  11  V,  vu(0)  >  10.6  V. 

Fig.  2  shows  the  tracking  performance  when  the  load  switches 


between  light  and  heavy,  with  yref 


1.5 

20 


.  Both  the  battery 


voltage  £  and  the  initial  supercapacitor  voltage  vu(0)  were  set  at  the 
nominal  value,  12  V  and  14  V,  respectively.  The  response  by  the 
averaged  model  is  plotted  with  red  curves  (in  the  web  version)  and 
the  response  by  the  SimPower  model  is  plotted  with  blue  curves. 
For  all  the  variables  v0,  ij,  and  12  by  the  SimPower  model,  the  curves 
are  saw-toothed  with  switching  frequency  20  kHz.  Thus  they  look 
like  solid  bands  and  the  width  of  the  band  shows  the  magnitude  of 
ripple.  For  the  output  voltage  v0,  the  band  is  thin  and  the  ripple  is 
small.  The  ripple  of  the  current  ii  is  visible  but  much  smaller  than 
the  ripple  of  12,  whose  peak-to-peak  value  is  nearly  5  A. 

During  the  first  0.2  s,  the  load  R  =  200  Q;  from  0.2  to  0.4  s, 
R  =  4  O;  from  0.4  to  0.6  s,  R  =  200  O.  The  battery  current  it,  is  similar 
to  i]  but  with  smaller  ripples  due  to  the  filter  capacitance  Q.  Also, 
the  supercapacitor  current  isc  is  similar  to  12  but  with  smaller  rip¬ 
ples.  During  light  load  period,  the  supercapacitor  absorbed  current, 
during  heavy  load  period,  it  supplied  more  than  5  A  of  current.  Both 
the  inductor  current  U  and  the  load  voltage  v0  were  tracked  under 
switching  load. 

Fig.  3  shows  the  response  to  load  switch  under  the  same  feed¬ 


back  gain  /<i  and  reference  yref 


1.5 

20 


but  £  is  reduced  to  11  V 


and  vu(0)  is  reduced  to  9  V.  The  reference  is  still  tracked.  The  most 
visible  difference  from  Fig.  3,  is  that  12  is  significantly  increased 


Response  under  K  ,  E=11V,  v  (0)=9V,  v  =20V,i  =1.5A 

r  1  u'  '  o.ref  b.ref 


time  (second) 

Fig.  3.  Response  to  load  switch:  £  =  11  V,  vu(0)  =  9  V,  feedback  gain  Ki. 


during  the  heavy  load  period.  The  steady  state  is  kept  at  1.5  A  as 
intended.  The  other  difference  is  that  the  overshoot  and  under¬ 
shoot  of  i]  and  v0  are  increased. 

Comparing  Figs.  3  and  2,  we  can  see  that  the  battery  current  and 
the  output  voltage  are  very  similar  under  different  state  of  charge  of 
the  power  sources.  The  supercapacitor  supplies  necessary  power  to 
protect  the  battery  and  to  maintain  the  desired  output  power. 

Although  robust  tracking  is  achieved  by  the  feedback  gain  I<h 
the  transient  dynamics  during  load  switch  is  not  satisfactory.  The 
overshoot  and  undershoot  of  v0  are  nearly  10%.  The  overshoot  and 
undershoot  of  ii  is  almost  100%. 


4.3.  Controller  for  smooth  transition  and  simulation 

To  reduce  the  overshoot  and  undershoot  of  both  ii  and  v0,  we 
designed  a  new  control  law  by  solving  the  optimization  problem 
(17).  By  choosing  77  =  50  and  w  =  1,  y  =  2000,  we  obtained  an 
optimized  feedback  gain  I< 2: 


-0.0056 

-0.0007 

-0.0001 

-0.0004 

-0.1011 

-0.0137 

-0.0014 

-0.0065 

0.0204 

-0.0650 

-7.8443 

1.0394 

-1.1275 

-7.6319 

We  used  this  feedback  gain  for  simulation  with  both  the  aver¬ 
aged  model  and  the  SimPower  model.  Different  values  for  the 
reference  yref  have  been  tested.  For  battery  voltage  and  super¬ 
capacitor  voltage  within  a  certain  range,  robust  tracking  have  also 
been  verified.  An  obvious  improvement  from  I< 1  is  that  both  the 
overshoot  and  undershoot  for  i'i,  v0  are  significantly  reduced. 

The  simulation  results  are  shown  in  Figs.  4  and  5. 

The  setup  for  generating  Fig.  4  was  the  same  as  that  for  gener¬ 
ating  Fig.  2,  except  that  the  feedback  gain  /<i  was  replaced  with  I< 2. 
As  compared  with  Fig.  2,  the  overshoot  and  undershoot  for  v0  was 
reduced  from  10%  to  less  than  5%.  The  overshoot/undershoot  for  ij 
were  almost  all  gone.  Similar  performance  improvement  can  be 
seen  in  Fig.  5. 


H.  Jung  et  al.  /  Journal  of  Power  Sources  267  (2014)  566-575 


573 


Response  under  E=12V,  v  (0)=14V,  vq  ref=20V,  ^  =1 .5A 


ldc-dc  con ve rte rs 


(Analog  interface 


Battery 


Micro¬ 

controller 


Supercapacitor 


Fig.  6.  Experiment  setup. 


Fig.  4.  Response  to  load  switch:  ii>ref  =  1.5  A,  v0,ref  =  20  V,  feedback  gain  K 2. 

5.  Experimental  results 

The  experimental  test  bench  was  constructed  at  a  reduced  scale, 
as  compared  to  some  real  applications.  The  circuit  is  described  in 
Fig.  1  and  the  parameters  are  given  in  Table  1.  For  the  battery 
module,  a  Panasonic  LC-R067R2P  was  selected.  The  supercapacitor 
consists  of  two  modules  of  BMOD0058  MAXWELL  58  F— 16.2  V  DC 
in  parallel.  Same  type  of  MOSFETs  were  used  for  the  bidirectional 
buck— boost  converters.  The  controller  was  implemented  with 
Texas  Instruments'  TMS320F28335.  The  switching  frequency  of  the 
dc— dc  converters  and  the  sampling  frequency  were  programmed  at 
18  kHz.  The  experiment  setup  is  shown  in  Fig.  6. 

The  state  feedback  control  u  =  I<2X  (I< 2  provided  in  (18))  was 
implemented  with  a  microprocessor.  The  initial  battery  voltage  was 
between  12  V  and  12.3  V,  and  the  initial  supercapacitor  voltage  was 
between  13.6  V  and  14.8  V.  The  load  was  switched  between  5  Q  and 
360  0. 

The  objective  is  to  track  the  load  voltage  v0  and  the  inductor 
current  i  1,  which  is  very  similar  to  the  battery  current  it,  except  with 


Response  under  K  E=12V,  v  (0)=14V,  vq  ref=16V^  ref=3A 


time  (second) 

Fig.  5.  Response  to  load  switch:  ilref  =  3  A,  voref  =  16  V,  feedback  gain  I(2. 


larger  ripples.  We  set  the  reference  for  v0  at  v0,ref  =  16  V,  18  V,  20  V, 
respectively.  For  the  reference  of  ip  we  tested  two  values 
ij  ref  =  1.5  A,  3  A.  To  verify  the  robust  tracking  performance,  same 
controller  11  =  K2X  was  used  for  all  tests. 

Fig.  7  is  the  experimental  response  to  the  nominal  reference 
ii.ref  =  3  A,  v0,ref  =  20  V.  Channel  1  (top,  blue  in  the  web  version, 
10  V/div)  is  the  load  voltage  v0.  Channel  2  (second  from  top,  red, 
2  A/div)  is  the  inductor  current  ij  for  the  battery  side  boost  con¬ 
verter.  Channel  3  (bottom,  green,  10  V/div)  is  the  supercapacitor 
voltage.  Channel  4  (third  from  top,  orange,  10  A/div)  is  the  inductor 
current  i2  for  the  supercapacitor  side  boost  converter.  The  super¬ 
capacitor  current  isc  is  very  similar  to  ij  but  with  smaller  ripples. 
The  time  scale  is  1  s/div.  Other  figures  have  the  same  channel  layout 
and  scales. 

For  the  initial  3.5  s,  the  load  was  360  Q.  Both  i  1  and  v0  tracked  the 
set  reference  values,  ij  ref  =  3  A,  vo  ref  =  20  V,  except  with  some 
ripples  and  noises.  At  about  3.5  s,  the  load  was  switched  to  5  Q.  At 
about  8  s,  the  load  was  switched  back  to  360  O. 

There  is  no  visible  change  in  ij  and  v0  after  each  load  switch.  In 
the  corresponding  simulation  results,  there  are  slight  overshoots/ 
undershoots  in  ij  and  v0  (less  than  1  V).  These  overshoots/un¬ 
dershoots  seem  to  be  dominated  by  the  ripples  and  noises  in  the 
experimental  results. 


Fig.  7.  Tracking  response  to  load  switch,  i,  ref  =  3  A,  v„  ref  =  20  V,  R  =  5  O,  360  Q. 


574 


H.  Jung  et  al.  /  Journal  of  Power  Sources  267  (2014)  566-575 


Fig.  8.  Tracking  response  to  load  switch,  i,  ref  -  1.5  A,  v0  ref  20  V,  R  5  fi,  360  Q 


The  load  switch  is  reflected  from  h  (Channel  4).  It  was  negative 
during  light  load  period  and  jumps  to  about  4  A  during  heavy  load 
period.  During  heavy  load  period,  the  supercapacitor  voltage 
slightly  decreases. 

This  result  verifies  the  robust  tracking  performance  to  load 
switch,  and  smooth  transition  for  battery  current  and  load  voltage, 
after  the  load  is  switched. 

To  verify  the  robust  tracking  performance  to  different  reference 
values,  we  obtained  the  results  shown  in  Figs.  8  and  9. 

Fig.  8  shows  the  tracking  response  to  ij  ref  =  1.5  A,  vo  ref  =  20  V. 
Both  outputs  were  robustly  tracked  with  smooth  transition  when 
the  load  was  switched  between  5  0  and  360  0.  The  difference  from 
Fig.  7  (h  ref  =  3  A)  is  that  the  supercapacitor  absorbed  less  current 
during  light  load  period  and  supplied  more  current  during  heavy 
load  period.  This  is  expected  since  the  role  of  the  supercapacitor  is 
to  supply  necessary  power  so  that  the  current  from  the  battery  is 
maintained  at  a  desired  value. 

The  decrease  of  the  supercapacitor  voltage  during  heavy  load 
period  is  more  visible  than  the  case  with  ij  ref  —  3  A.  Because  of  this, 
the  supercapacitor  current  slightly  increased  as  it  was  discharged. 
Another  difference  from  Fig.  7  is  that  during  heavy  load  period,  the 


Fig.  9.  Tracking  response  to  load  switch,  ilref  =  3  A,  voref  =  16  V,  R  =  5  Q,  360  Cl. 


ripples  and  noises  of  h  and  v0  were  larger.  The  supercapacitor 
voltage  also  had  larger  ripples. 

Fig.  9  shows  the  tracking  of  h.ref  =  3  A,  vo  ref=  16  V.  As  compared 
with  Fig.  7,  the  output  power  demand  is  smaller  so  the  power 
supplied  by  the  supercapacitor  is  reduced,  as  reflected  by  the 
smaller  z'2. 

We  tested  other  combinations  of  output  reference,  battery 
voltage  and  supercapacitor  voltage  and  verified  that  reference 
tracking  can  be  achieved  within  a  large  range  of  parameters.  The 
tracking  responses  were  very  similar.  The  output  voltage  and  bat¬ 
tery  current  were  always  around  the  desired  reference  values.  All 
the  changes  of  reference  values,  the  source  voltages  and  load 
resistance  were  handled  by  the  supercapacitor. 

6.  Conclusions 

We  addressed  some  control  design  problems  for  a  power  system 
driven  by  battery/supercapacitor  hybrid  energy  storage  devices. 
The  control  objectives  include  tracking  the  references  for  two  var¬ 
iables  in  the  presence  of  uncertain  conditions  and  achieving 
smooth  transition  during  load  switch.  Based  on  a  newly  derived 
state-space  averaged  model,  the  control  design  problems  were 
converted  into  optimization  problems  with  LMI  constraints.  The 
effectiveness  of  the  control  design  methods  was  validated  by  both 
simulation  and  experiment.  The  results  developed  in  this  paper  will 
find  possible  applications  in  electric,  hybrid  and  plug-in  hybrid 
electric  vehicles,  wind  systems,  photovoltaic  systems,  and 
microgrids. 
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